Kohn Anomalies and Electron-Phonon Interaction in Graphite 
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We demonstrate that graphite phonon dispersions have two Kohn anomalies at the r-E2 9 and 
K-Al modes. The anomalies are revealed by two sharp kinks. By an exact analytic derivation, we 
show that the slope of these kinks is proportional to the square of the electron-phonon coupling 
(EPC). Thus, we can directly measure the EPC from the experimental dispersions. The r-E2 9 and 
K-A'j EPCs are particularly large, whilst they are negligible for all the other modes at T and K. 

PACS numbers: 63.20.Dj, 63.20.Kr, 71. 15. Mb, 78.30.-j 
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Carbon nanotubes (CNTs) are at the core of nanotech- 
nology research. They are prototype one dimensional 
conductors. Metallic nanotubes are predicted to be one- 
dimensional quantum wires with ballistic electron trans- 
port 0. However, high field electrical transport mea- 
surements show that the electron-phonon scattering by 
optical phonons at K and T breaks down the ballistic 
behavior |2|. Electron phonon coupling (EPC) is thus 
the fundamental bottleneck for ballistic transport. Ra- 
man spectroscopy is a most used characterization tech- 
nique to identify CNTs in terms of their size and elec- 
tronic properties [3j . The optical phonons at K and T are 
the phonons responsible for the Raman D and G peaks 
in carbons Q. The frequency and the intensity of the 
Raman modes are determined by the EPC matrix ele- 
ments |5J. The determination of the EPCs is necessary 
to settle the 35 years debate on the nature of the Raman 

D peak in carbons HSUHHISHIUEl- Fmailv ; 
although graphite phonon dispersions have been widely 
studied, several contras ting th eoretical dispersions were 
proposed [fj IjJ, la |9j, llIJ, LLii LLJ] . In particular, the origin 
of the large overbending of the K-A^ branch is not yet 
understood and is associated to an intense EPC j lL Ivjl. 1 1 ( )l] . 
In principle, the electronic and vibrational properties of 
CNTs can be described by folding the electronic and 
phonon dispersions of graphite. The precise determina- 
tion of the graphite EPCs is thus the crucial step to un- 
derstand the properties of any carbon based material and 
CNTs in particular. It is then surprising that, despite the 
vast literature on carbon materials and CNTs, no experi- 
mental determination or first principle calculations of the 
graphite EPCs has been done so far, to the best of our 
knowledge. 

Here we show that in graphite the EPC matrix ele- 
ments at r and K can be directly extracted from the 
phonon dispersions. We demonstrate two remarkable 
Kohn anomalies in the phonon dispersions at T and K, by 
an exact analytic derivation and accurate density func- 
tional theory (DFT) calculations. We prove that the 
slope of the anomalies is proportional to the EPC square. 

A key feature of graphite is the semi-metallic character 
of the electronic structure. In general, the atomic vibra- 



tions are partially screened by electrons. In a metal this 
screening can change rapidly for vibrations associated to 
certain q points of the Brillouin Zone (BZ), entirely de- 
termined by the shape of the Fermi surface. The con- 
sequent anomalous behavior of the phonon dispersion is 
called Kohn anomaly Kohn anomalies may occur 

only for wavevectors q such that there are two electronic 
states ki and k.2 = ki + q both on the Fermi surface . 
The electronic bands dispersions of graphite are, essen- 
tially, described by those of an isolated graphene sheet. 
In graphene, the gap between occupied and empty elec- 
tronic states is zero only at the two equivalent BZ points 
K and K'. Since K' = 2K, the two equivalent K points 
are connected by the vector K. Thus, Kohn anomalies 
can occur for q = T or q = K. 

We perform calculations within the generalized gradi- 
ent approximation (GGA) using the density func- 
tional perturbation theory (DFPT) scheme 0]> which 
allows the exact (within DFT) computation of phonon 
frequencies at any BZ point. We use the plane- waves 
(90 Ry cut-off) and pseudopotential [18j approach. The 
electronic levels are oc cup ied with a finite fictitious elec- 
tronic temperature a |l9j |. This smears out the discon- 
tinuities present in the Fermi distribution for a = 0, and 
the exact result is recovered for a—tO. The experimental 
lattice (a oxp = 2.46 A, c = 6.708 A) is used for graphite, 
while for graphene we consider both the graphite a OX p 
and the theoretical values (a t h = 2.479 A). Graphene 
layers are separated by 7.4 A of vacuum. 

Fig. compares the measured optical branches with 
our calculations at a = 0.02 Ry. Phonon frequencies 
are computed exactly for a series of points along the 
high symmetry lines jT-K and jT-M and then interpo- 
lated with a spline. The agreement with experiments is 
~ 2%, which is the expected accuracy of DFT-GGA. At 
r the experimental dispersion is very well reproduced by 
the calculations with a cxp . At K the upper branch is 
better described by the calculations with ath- 

The most striking feature of these dispersions is the 
discontinuity in the frequency derivative of the highest 
optical branches (HOB) at T and at K (F,2 g and A[ 
modes). Indeed, near T, hujq = org + fkov + 0(q 2 ), u>q 
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FIG. 1: (Color online) Upper panel: Lines are the phonon 
dispersion of graphene (GE), calculated at the experimental 
and equilibrium lattice spacings (a exp and ath)- Experimen- 
tal data from Ref. 0. The red straight lines at T and K are 
obtained from Eqs. 1811(71 The two lower panels correspond to 
the dotted window in the upper panel. Here, graphite (GI) 
computed frequencies are also shown. The points are theoret- 
ical frequencies obtained by direct calculation. A single GE 
band corresponds to two almost-degenerate Gl-bands. 



being the phonon frequency of the HOB at the q wavevec- 
tor. Similarly, near K, ?iwK+q' = ok?' + + 0{q' 2 ). 
Such dependencies cannot be described by a finite set of 
interatomic force constants, or by a set decaying expo- 
nentially with the real-space distance. In these cases the 
dynamical matrix dependence on q would be analytic, 
and, because of symmetry, the highest optical branch near 
r and K would have a flat slope («r = aK = 0). Thus, a 
non zero ctr or ajc indicates a non-analytic behavior of 
the phonon dispersion, due to a polynomial decay of the 
force constants in real space. This explains why it is im- 
possible for any of the often used few-nearest-neighbors 
force constants approaches to properly describe the HOB 
phonons near K and T H E H El H Thc 
graphite HOB are almost indistinguishable from those 
of graphene, Fig.^ The non-analytic behavior at T and 
K is also present in graphite. At T the HOB is doubly 
degenerate, consisting of in-plane anti-phase E2 g move- 
ments. Near T the two modes split in a upper longitu- 
dinal optical (LO) branch and a lower transverse optical 
(TO) branch. «r only for the LO branch. 

Fig. plots the HOB as a function of a, to clarify the 
nature of the discontinuities. The results for ct = 0.01 Ry 
and <t = 0.02 Ry are similar, indicating that, on the scale 
of the figure, the a — > limit is reached. Increasing a, 
the non-analytic behavior is smoothed out. Particularly 
striking is the behavior around K, where for a > 0.20 Ry 
the dispersion is almost flat. Within DFPT, the smearing 
a affects virtual transitions between occupied and empty 
states, differing in energy by < a. Thus, the smooth- 
ing of Fig. [21 indicates that the HOB discontinuities are 
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FIG. 2: (Color online) Graphene HOB around F and K as a 
function of smearing. Points are calculated frequencies for ath • 
Lines guides to the eye. The red straight lines plot Eqs. I8I1UI 



Kohn anomalies 0] , since they are related to an anoma- 
lous screening of the electrons around the Fermi energy. 

We compute the EPC matrix elements, to understand 
why the Kohn anomalies affect only the HOB and not 
the others. For a given phonon mode at the reciprocal- 
space point q, we call AVq and An q the derivative of 
the Kohn-Sham potential and of the charge density, with 
respect to a displacement along the normal coordinate of 
the phonon. We define: 

(k + q,t|AV q [AnJ|kjWty(2Mo; q ), (1) 

where we consider explicitly the dependence of AVq 
on An q , as indicated by AVq[An q ]. |k, i) is the elec- 
tronic Bloch eigen-state with wavevector k, band in- 
dex i, and eigen-energy ek,i- M is the atomic mass, 
and g is an energy. The dimensionless EPC 
is A q = 2J2k,i,j l5(k+ q ) i ,k J | 2 <5(ek+q, l - £F)<5(ek,j - 
€F)/(fujj< 1 NFNk), where J^k 1S a sum on BZ vectors, 
Nf is the density of states per spin at the Fermi energy 
€f- In graphene the Fermi surface is a point, JVf is zero 
and A n is not well-defined. Thus, we evaluate 



2(5, 



VVp 

Jo ' 



where (■ ■ ■ )f indicates the average on the Fermi surface of 
l.9(k+ q ) 4 ,kj| 2 , and J q = 1/N k T,i,j,k 5 0k+q,i - e F )S{e k j - 
e F ). In graphene, (g^) F = Y,i,j\9(2K)i,Kj\ 2 /4, and 
(9f)f = Y^lj ISKi.Kj | 2 /4, where the sums are performed 
on the two degenerate it bands at eF- For the HOB, we 
obtain (g r ) F = 0.0405 eV 2 [HI, and (g^} F = 0.0994 eV 2 , 
corresponding to 2(5 2 )p/(fia; q ) of 0.41 eV and 1.23 eV at 
r and K, respectively. 2(<7 2 c )p / (Swk) is much smaller 
(0.02 eV) for the doubly degenerate 1200 cm -1 phonon 
at K and zero for all the other phonons at K and at 
r, consistent with the absence of Kohn anomalies for 
all these branches. The EPCs values for the graphene 
HOB are very large in absolute terms. They are com- 
parable to the 39 K superconductor MgB2, for which 
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AqiVp/Jq = 1.6 eV [22| for each of the two strongly- 
coupled E2 g branches at A. On the contrary, in alka- 
line doped C 6 o, the largest 2{g 2 l ) F /(fkj q ) for individual 
phonon modes are much smaller (0.02-0.034 eV |2(J). 

We consider the expression of w q according to pertur- 
bation theory (within DFT), to understand the ab- 
sence of a Kohn anomaly for the TO E2 g phonon at T and 
to correlate the constants or and «k with the magnitude 
of the EPC. ujq = -J 1 Dq/M, where D q is the dynamical 
matrix projected on the phonon normal coordinate: 



4 y, Kk + q, e |Al/ q [An q ]|k, )| 2 



iVk 



J Ari;(r)if(r,r')An q (r') d\ dV 



+ / n(r)A 2 V b (r) d\. 



(2) 



Here a factor 2 accounts for spin degeneracy, ^ c is a 
sum on occupied and empty states, n(r) is the charge 
density, K(r,r') = S 2 E Hxc [n]/(Sn{r)5n(r')), E Hxc [n] 
is the Hartree and exchange-correlation functional, and 
A 2 V b is the second derivative of the bare (purely ionic) 
potential. From previous considerations, 
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If the dependence of £> q on q were analytic over all the 
BZ, a r = a K = (E.g. D q - D T = 0{q 2 )). For q^T 
and K, the energy denominators in the sum of Eq.[3]arc 
finite and D q is analytic. The energy denominators go 
to zero for q = T (when k = K or k=2K) and for q = K 
(when k = K), when o and e correspond to the it and tt* 
bands near the Fermi energy. Due to these singularities, 
the dynamical matrix is non- analytic for q = K, T, thus 
ar and ckk can be different from zero. To compute a q , 
we can replace in Eas. l3l4l the full dynamical matrix Z? q 
with its non-analytic component D q , which includes only 
the sum of Eq. [21 restricted to the it electronic bands in 
an arbitrarily-small but finite circle of radius k m around 
the Fermi surface at the K and/or 2K points 23]. 
For q near T, using the definition of Eq. Q 
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lot 



k'<k„ 



, 2 ,/ lff(K+k'+q)7r«,(K+k')7iT 
£K+k',7r — £K+k'+q,7r* 



(5) 



Here we have used the substitution l/A^^k = 
V3/2 J d 2 k' (the graphene BZ area is 2/V3), q and k 
points are in units of 2ir/a, a being the lattice spacing, 
k' = k K. A factor 2 accounts for the contribution of 
the two equivalent Fermi points. For small q and k', the 
EPC matrix elements in the numerator in Eq. 0arc: 

LO 00 

l.9(K +k ' + q)^,(K + k J = <3r>F [1 ± co S (9 + 9% (6) 



where is the angle between k' and q, and 9' is the angle 
between k' + q and q [24| . The + or — sign refers to the 
LO and the TO modes, respectively. In graphene, the 
electronic bands near the Fermi level have a conic shape. 
Therefore, for small q and k', 



£K+k',7r — CK+k'+q.Tr* 



-13k' - 0\\s! - 



(7) 



where = 14.1 eV is the slope of the tt bands within 
GGA. Replacing Ea . 15161 and in Eq. [31 one obtains: 



LO 
, TO 



lim - 

q^O q 



4v / 3( 3r >F 




d 2 k' 



l±cos20 l±cos(0 + 0') 



2k' 



k'<k n 



4V3(, 9r ) F 





dy Jd0 





fc' + |k' + q| 
y ± ycos(0 



y - 



y 



2y cos 



where y = k/q, and 0' = arctan[y sin 0/(1 + ycos0)]. 
The integral is zero for the TO mode and ir 2 /4 for the 
LO mode. Thus, as expected, a™ = and 







V3tt 2 = 397 cm -1 . 



(8) 



For the K point, we use an analogous procedure, for tran- 
sitions from points in the neighborhood of K to points in 
the neighborhood of 2K. The EPC matrix elements are: 



l.9(2K+k'+q')*-,(K+k>| 2 = (.9 K >F (1 +COS0"), (9) 



with 9" the angle between k' and k' + q' 
becomes: 



2V3(. 9 2 C ) F 1 
a K = ^ — nm — 

q'^o q' 



d\' 



1 

fc 7 



1 



I2J- Eq.H 

COS0" 



fe' + lk' + q'l 



fe'<fc m 

The limit of the integral is 7r 2 /2, and therefore 







'.y/STT 2 = 973 cm" 1 . 



(10) 



The resulting linear dispersions are plotted in Fig.^ an d 
Fig. [21 As expected, the phonon slopes near the discon- 
tinuities are very well reproduced. Finally, we note that, 
within a first-neighbors tight-binding (TB) approxima- 
tion, the EPCs at T and K arc not independent, since: 



(■9k) F 

(ff r >F w r 



(11) 



The validity of this relation is supported by our DFT re- 
sult, for which ((gj^F WK)/((g r )p w r ) = 2.02. 

Eas. 181101 allow us to directly measure the EPCs at T 
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and K from the experimental phonon slopes. The phonon 
branches around T have been determined by several 
groups with a close agreement of the measured data |13L 
The fit for q < 0.15 of the measurements in Fig. ^ 9] 



with fku„ = fuvr 



*r°<L+ "f1 2 y S ives a r° = 340 cm -1 , 
and (ffr)F = 0.035 eV 2 |2l|. The available data around 
K are much more scattered. However, from Eq. 1111 we 
get (#k)f = 0-086 eV 2 . These values are in excellent 
agreement with our calculations and validate our results. 

The EPCs near K between it* bands 25] allow the ac- 
curate determination of the intensity and shape of the 
Raman D peak. This will be the topic of future pub- 
lications. Here, we remark that the branch has the 
biggest (5k)f amongst K phonons. Thus, the D peak is 
due to the highest optical branch starting from the K-A^ 
mode |j, |9|, llOl I 111, not to the branch starting from the 
K-E' mode, as m|E IS 0, 11 E3 • Also > the D P eak shifts 
linearly with laser excitation energy (^50 cm _1 /eV [l2jh 
This is at odds with the flat slope of the K-A^ branch 
obtained by previous calculations 0, 0, S EI EH • 
But it is consistent with the K-A'j linear-dispersion we 
get. The D peak dispersion reflects the slope of the Kohn 
anomaly at K and provides another independent mea- 
sure of the EPCs. From ref. [5j, for q' — > 0, the D peak 
dispersion is ~ ok//3. This gives (<?k)f ~ 0.072 eV 2 
and, from Ea.HT| (pf,) F - 0.029 eV 2 . The EPCs derived 
in this way are a lower limit since the experimental D 
peak dispersion is measured with laser excitation ener- 
gies > 1 eV 0. This corresponds to the phonon slope 
at q' > 0.035, which, from Fig.0 underestimates by 30% 
the slope at q' = 0. Taking this into account, the EPCs 
independently inferred from the D peak dispersion are as 
well in excellent agreement with our calculations. 

Due to the reduced dimensionality, we predict even 
stronger Kohn anomalies for metallic CNTs, and no 
anomaly for semiconducting CNTs. This is the key to dif- 
ferentiate the electrical nature of CNTs by Raman spec- 
troscopy. A softening of CNT phonons corresponding to 
the graphene T-E 2g mode was recently reported 0] . We 
expect a stronger softening for the phonons correspond- 
ing to the graphene K-A^ mode, since (gK)F>(3r)F- 

In conclusion, we have demonstrated the presence of 
two remarkable Kohn anomalies in the phonon disper- 
sions of graphite, revealed by two kinks for the r-E2 9 
and K-A'j modes. Even if Kohn anomalies have been 
observed in many materials [TtI ]. graphite is the first 
real material where a simple analytic description of the 
anomaly is possible. Indeed, we proved, by an exact an- 
alytic derivation, that the slope of the kinks is propor- 
tional to the ratio between square EPC matrix elements 
and the slope of the tt bands at the Fermi energy. As 
a consequence, we directly derived the EPCs at T and 
K from the experimental phonon dispersions. The val- 
ues we obtain are in excellent agreement with calcula- 
tions. Finally, our (<?p)f and (<7 2 c )f values, with Eas. 16191 
and [25j, can be used to determine the mean free path 



for electron scattering by optical phonons. This gives the 
limit of ballistic transport in CNTs and is of great sci- 
entific and technologic importance Q. This calculation 
can be done, within the folding model, using the Fermi 
golden rule, and will be reported elsewhere. 
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